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Abstract 

The gravitational field of a black hole is strongly localized near its horizon when the 
number of dimensions D is very large. In this limit, we can effectively replace the black 
hole with a surface in a background geometry {e.g., Minkowski or Anti-deSitter space). 
The Einstein equations determine the effective equations that this ‘black hole surface’ (or 
membrane) must satisfy. We obtain them up to next-to-leading order in 1/D for static 
black holes of the Einstein-(A)dS theory. To leading order, and also to next order in 
Minkowski backgrounds, the equations of the effective theory are the same as soap-film 
equations, possibly up to a redshift factor. In particular, the Schwarzschild black hole is 
recovered as a spherical soap bubble. Less trivially, we find solutions for ‘black droplets’, 
i.e., black holes localized at the boundary of AdS, and for non-uniform black strings. 
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1 Introduction 


Recently it has been demonstrated that black hole physics can be efficiently solved in an 
expansion around the limit of large number of dimensions, D ^ oo [1, 2, 3, 4, 5, 6 , 7]. In 
this limit there appears a near-horizon region that is universal for all non-extremal neutral 
black holes, and which encompasses the small radial extent ~ tq/D of the gravitational 
field outside of a horizon of radius ~ rg [3] . One important feature is that the quasinormal 
spectrum of the black hole splits into modes of high frequency, w ~ D/vq, and of low 
frequency co ~ 1/ro [4, 6 , 7, 8 ]. The latter are particularly interesting, since they are 
fully localized in the near-horizon region (where they are normalizable excitations) and 
are decoupled from the asymptotic ‘far zone’ to all perturbative orders in 1/D. This split 
into two scales makes it natural to try to integrate out the high-frequency, short-distance 
degrees of freedom to obtain a fully non-linear effective theory of the long-wavelength, 
decoupled dynamics. Such a theory should capture the physics of black holes on lengths 
and timescales S> tq/D, hence allowing for fluctuations on scales comparable to the horizon 
radius, ~ tq, or (as appropriate for black branes [ 1 , 2 ]) rsj ro/^/D. 

We can also motivate this effective theory in a more geometric fashion. In the limit 
D —>■ oo the gravitational field of the black hole vanishes outside the near-horizon region, 
and thus there is a neat separation between the black hole and the background where 
it resides (e. 5 ., Minkowski or (Anti-)deSitter spacetime). The black hole can then be 
effectively identified with a particular surface in this background. Since ultimately the 
properties of the black hole are dictated by Einstein’s equations, it must be possible to 
derive from them a set of equations that must satisfy. These constitute the effective 
theory of the black hole in the large D limit. 

The black hole is then described by a set of ‘collective coordinates’ which specify 
the embedding of in the background spacetime, and which vary over scales much 
larger than tq/D, where tq is a characteristic length of The approach to obtain the 
effective theory employs the parametric separation between the large radial gradients and 
the smaller temporal and spatial gradients along the horizon, which allows to solve the 
radial dependence of the Einstein equations. Then, the vector-constraints in the radial 
direction yield the effective equations for the embedding functions. Readers familiar with 
other effective theories of black holes, in particular the fluid/gravity correspondence [9] 
and the blackfold approach [10], will recognize similarities here. They are all based on 
a parametric separation between the fluctuations that are transverse to the horizon and 
those that are parallel to it. However, since D is a parameter of the theory instead of a 
parameter specific to some solutions, in principle the 1/D expansion allows to tackle a 
larger set of problems —at the expense of possibly losing accuracy at finite values of D or 
missing phenomena which are non-perturbative in 1/D. 

This effective theory is an important step in the program to understand gravity in 
the large D limit. General Relativity in vacuum (possibly with a cosmological constant. 
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and without compactified dimensions) is essentially a theory of black holes and gravita¬ 
tional waves. The large D effective theory is a reformulation of the black hole sector of 
General Relativity in terms of membrane-like variables. The decoupling property of the 
near-horizon zone implies that, as long as its gradients remain much smaller than Djro, 
the effective membrane can not radiate any gravitational waves to the far zone at any 
perturbative order in the expansion. Conversely, no gravitational waves from the far zone 
of frequency <C D/tq can interact with the effective membrane; and, while far-zone waves 
of frequency ~ D/tq or larger can penetrate to the near zone, they are perfectly absorbed 
by the horizon on a short timescale ~ tq/D and do not interact with the low-frequency 
modes of the effective membrane [2]. So the two sectors of the theory —black holes and 
gravitational waves— appear to be well separated, with the coupling between them being 
non-perturbative in 1/D. However, there do exist black holes that when Z) —)• oo have 
large spatial gradients ^ D/vq along their horizons. Moreover, in the evolution of some 
horizons it can occur that initially small temporal gradients grow to values ~ D/vq. Such 
situations imply breakdowns of the applicability of the effective theory, and are reminiscent 
of the breakdown of hydrodynamics when turbulence develops. 

In this article we start to develop the large D effective theory of black holes by fo¬ 
cusing on the simplest case of static configurations of neutral black holes, possibly with a 
cosmological constant. The equations for the embedding of in the background take a 
remarkably simple form: if the trace of the extrinsic curvature of is K, and the redshift 
factor on Ss is gu, then the effective equation at leading order is 


V-9tt K = 2k, (1.1) 

where the constant n gives the surface gravity of the black hole. When this is satisfied, our 
results provide explicitly the near-horizon black hole metric that solves the full Einstein- 
(A)dS equations in the leading large D limit. Observe that in backgrounds like Minkowski 
space where there is no redshift, this is the same as the equation for a soap film. The 
stress-energy tensor of the effective membrane that lives on T,b is also simple: it is a 
modulation along the membrane directions of the quasilocal stress tensor of the large-T) 
black hole. 

We also obtain the next-to-leading order corrections to the effective theory. Although 
in general the form of the equations becomes a little more complicated than (1.1), in 
the absence of redshifts we get the same equation K = const, at the next order in the 
expansion. 

In order to test the usefulness of the theory, we have applied it to obtain several 
non-trivial solutions. Staticity greatly restricts the number of possible black holes. For 
instance, in all dimensions the unique asymptotically flat, static black hole of Einstein’s 
vacuum equations is the Schwarzschild-Tangherlini solution [11], which is recovered from 
(1.1) as a spherical soap bubble — correspondingly, the spherical soap bubble is the unique 
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closed surface of constant mean curvature in flat space in any dimension^. However, we 
will show that the effective theory easily yields other static solutions: black holes localized 
at the boundary of AdS (‘black droplets’) [13], and non-uniform black strings in asymp¬ 
totically flat space [14, 15, 16]. The construction of these solutions at specific finite values 
of D required sophisticated numerical solution of systems of partial differential equations. 
In contrast, the large-H equations for these problems are single ordinary differential equa¬ 
tions that can be easily solved, when not in an analytical exact or approximate form, at 
least with a one-line command of NDSolve in Mathematica. 

In the next section we introduce a formalism that is adequate for the resolution of the 
problem. Then in sec. 3 we solve the Einstein equations at leading order in 1 /D and derive 
the effective equation for the surface We also discuss simple examples, and obtain 
the stress-energy tensor for the effective membrane on In sec. 4 the effective equation 
is solved to obtain black droplets in AdS. Sec. 5 contains the derivation of the effective 
theory at next-to-leading (NLO) order. This is then applied in sec. 6 to the construction 
of non-uniform black strings. Sec. 7 finishes the paper with some brief remarks. The 
appendices contain technical details and elaborations on asides mentioned in the main 
text. 


2 Set up 


In our metric ansatz we separate a radial direction p, where gradients are of order H, from 
all other spacetime directions along which variations are smaller. To this effect, we use a 
decomposition in ‘synchronous gauge’ 

ds'^ = + gf,^{p,x)dx^dx^ , ( 2 . 1 ) 


where p, u run over all the D — 1 directions orthogonal to p. The Einstein equations in 
vacuum, with a cosmological constant 


{D-l){D-2) 

2£2 


( 2 . 2 ) 


can be written in terms of the intrinsic, and extrinsic, K^y, curvature tensors of the 
{D — l)-dimensional constant-p surfaces, 


- K^^yK\ = R + 


{D-l)[D-2) 

£2 


v.a:;: - = 0 

= Rf^y + 

A Jy 

K^u = ^g^^d,g^y. 


'^This is widely believed among differential geometers to be the case, see e.g., [12]. 
Gibbons for this reference. 


(2.3) 

(2.4) 

(2.5) 

( 2 . 6 ) 

We thank Gary 
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Eqs. (2.3) and (2.4) are respectively the scalar and vector constraints, while (2.5) is the 
‘dynamical evolution’ equation in the radial direction. Knowing eq. (2.6) can be 

integrated to obtain the metric. It is convenient to consider separately the equation for 
K obtained from eq. (2.5), 

+ K^ = (D- lflC\p, x) , (2.7) 

where we define IC{p, x) by 

■ P-®' 

We take the metric to be static. In the spatial directions orthogonal to p, we have non¬ 
trivial dependence on a number p = 0{D^) of spatial directions z°’, which are orthogonal 
to a n -b 1-dimensional space that we take to be a sphere We introduce 

n = D-p-3, (2.9) 


which can be used as the large expansion parameter instead of D. Our metric ansatz is 
then 


ds^ = N‘^{p, z)- 


dp"^ 


{D - 1)2 ~ ^ Sab{p, z)dz°‘dz^ + R2(p, z)qijdx"dx^ (2.10) 

where a = 1 ,... ,p, and qij is the metric on the unit 

The different metric functions will be assumed to scale with n in specific ways. In 
order to get oriented, note that one solution that we intend to recover is the leading-order 
near-horizon geometry of the Schwarzschild black hole [3], which can be written in the 
ansatz (2.10) with p = \ as^ 

ds^ = —Tq tanh^(p/2)(it^ -b rg -b — lncosh(/3/2)^ -b sin^ z (in„+i)(2.11) 

Observe that in the sphere radius we are keeping terms of order 1/n: due to the large di¬ 
mensionality of the sphere, such terms enter (through traces) in the leading order equations 
and thus must be kept at this order. The takeaway here is that, while the p-dependence 
in V appears at the leading order, instead in gab and R it is at C>(l/n). 

We then assume that 


K\,K = 0{n), K\,K^, = 0{1), R = 0{n^). 

We also assume that N, V, R are 0(1) and that 

gab, SaR = 0(1), or 0(l/n). 

In these two cases the Ricci tensor of gab is, respectively, 

= 0(1), or 0(n). 


( 2 . 12 ) 

(2.13) 

(2.14) 


The first instance in (2.13) will be exemplified in sec. 4, and the second one in sec. 6. 

^This can be extended to other spaces with intrinsic curvature oc at large n. 

®Here p is twice the one in [3]. 
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3 Effective theory: leading order 

Our strategy is to first solve eqs. (2.3) and (2.5) for the radial dependence of the extrinsic 
curvature with regularity at the black hole horizon, and then obtain the metric from (2.6). 
Radial integrations leave an undetermined function of 2 ;“ which cannot be eliminated by 
gauge choices: this is the one collective degree of freedom of the static black hole. The 
vector constraint (2.4) then yields a non-linear differential equation, on 2 ; only, for this 
degree of freedom, which is the effective equation we seek. 


3.1 Solving the leading order eqnations 

According to our assumptions we can separate the leading-order, p-independent terms as 


N{p,z) 

= No{z) + 0(l/n), 

(3.1) 

gabiP, Z) 

= labiz) + 0{l/n), 

(3.2) 

R{p, z) 

= TZiz) + 0{l/n), 

(3.3) 

IC{p, z) 

= ^ + 0(l/n). 

ro{z) 

(3.4) 


In this case we can integrate (2.7) to find^ 

The divergence at p = po{z) is the expected pole at the horizon, coming from K^t (all other 
components of the extrinsic curvature must be regular there). The metric to leading order 
is invariant under p —)• p + f{z), which we use to set po{z) = 0. Similarly, by rescaling p 
by an 0(1) function of 2 ;“ we can reach a gauge in which, to leading order, 

Nq{z) = ro(2;). (3.6) 


Then 


n 


K = —^ coth p. 
royz) 


Next, the equation for K^t is 


n 


ro{z) 


dpK\ + KK\ = 0, 


since BJ't + can be neglected as it is of lower order. This equation is solved as 

n 


K\ = 


(3.7) 


(3.8) 


(3.9) 


ro( 2 ;) sinhp’ 

where the integration function of 2 ; has been fixed again by requiring that K^t have a pole 
at p = 0. 


'^Henceforth in this section we omit the symbol 0{l/n) to unburden the notation. In addition, the 
subsequent analysis is valid only when rp > 0. 
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Since and K'^j are 0(1), in order to solve for them we need R°^b and R^j. The 
components of the curvature in a constant-p section 

g^udx^dx^ = —\/^{p, z)dt^ + ^ab{z)dz'^dz^ + TZ^{z)qijdx^dx^ (3.10) 


are readily obtained (see appendix A). The scalar curvature gets negligible contributions 
from the time direction, and is given at leading order by 


7^ = n2^(l-(D7^)2) , 

(3.11) 

where we abbreviate 


(D7^)2 = 7“^aa7^5b7^. 

(3.12) 

Comparing (3.11) to eqs. (2.8) and (3.6) we find 


rl) = P + wdl - mf) ■ 

(3.13) 


It is easy to see that this equation is actually equivalent to the scalar constraint (2.3) at 
leading order. 

Since the scalar curvature is dominated by the components along the sphere we have 


Using the covariant derivative for the metric 7 ^ 6 , the curvature along the direc¬ 
tions, obtained from (A.3), is 

. (3.15) 

Here we take into account that in the second case in (2.13), (2.14) the intrinsic curvature 
of 7 af, can contribute at the same order. 

Now we have all the terms needed to integrate the equations from (2.5), 


71 71 

-dpK\ + KK\ = R\ + , 

ro 

71 • ■ 71 ■ 

-dpK^^+KK^,=R}, + -^5^,. 

Imposing regularity at p = 0 we find 

K'^b = ?’o/“6tanh(p/2), 

K^j = — tanh(p/2), 

ro 

where we have defined the tensor 

f f X lab{z) DaD^,7^ ^'^'>Rab 

fabiz) = - — + — 


(3.16) 

(3.17) 

(3.18) 

(3.19) 


(3.20) 
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It is now straightforward to integrate the extrinsic curvatures (3.9), (3.18) and (3.19) 
to obtain the metric components. This gives 

V(p, 2 ) = Vo(2;)tanh(p/2), (3.21) 

gab{p,z) = Jabiz) + -rl{z) fabiz) lncosh(/9/2), (3.22) 

n 

R{p,z) = TZ{z) ^1 H—lncosh(/9/2)^ . (3.23) 

Vo{z) is an integration function from the p-integration of K^f For the ab and ij components 
the integration functions of z have been absorbed in a C>(l/n) redefinition of ^ab{z) and 
TZ{z). However, Vo{z) cannot be absorbed in that way. 

Up to this point we have solved all the radial dependence of the metric and it only 
remains to impose the vector constraint (2.4), whose only non-trivial component is along 
z“. To leading order it takes the form 

TitK^ + + (^fe In ^) K\ + = 0 . (3.24) 


Plugging in our previous results and using the identity 

n^n 7? \ 1 

(dhimZ) =^(41n7^-41nro) , (3.25) 

derived from (3.13), we find, 


dalnVoiy) - 5alnro(z) 



n 


(dblnTZ) sinh^(/9/2) = 0 . 


(3.26) 


Under either of the two cases in (2.13), (2.14), the last term is subleading and can be 
neglected. Thus, consistently, the p-dependence cancels out of the leading-order equation, 
which requires that Vq be proportional to tq. Observe that this condition is equivalent to 
requiring that the surface gravity k be uniform on the horizon. Indeed, the precise relation 
is 


Vo{z) = — ro(z). 
n 


(3.27) 


Defining a rescaled surface gravity 


(3.28) 


n 

our solution of the Einstein equations, valid in the near-horizon region to leading order in 
1 /n, is 


ds^ 


rQ(z) 4k^ tanh^(/3/2)(it^ -|- 

\ n‘‘‘ 

+V?‘{z) (cosh(/9/2))" d^n+i ■ 


n 


+ ( 7abiz) + -ro(z)/afe(2)lncosh(/9/2) ) dz°‘dz^ 


(3.29) 


Although K could be absorbed by rescaling the time coordinate, this may not be convenient, 
since the normalization of t is fixed by matching to the far-zone. 
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The geometry (3.29) can be interpreted as a modulation along of a near-horizon 
Schwarzschild black hole solution. The functions ro(^), 7ab(-z) and TZ{z) vary slowly com¬ 
pared to radial gradients, and one of them, say 'R-{z)^ can be regarded as the single 
collective degree of freedom for the black hole. Indeed, it is possible to give an alternative 
derivation of the effective theory following these ideas, as described in appendix B. 


3.2 Effective equation 

The metric (3.29) in the near-horizon geometry must be matched to the far-zone back¬ 
ground in the common ‘overlap zone’, which corresponds to 1 <C p <C n in (3.29). The 
matching must be such that the metric induced on a constant-p surface there, namely, 

—VQ{z)dt^ + 'Jabiz) dz°'dz^ + TZ‘^{z)dfln+i, (3.30) 

and its extrinsic curvature, are the same when approached from either zone. From (3.21) 
and (3.7) we see that 

^\eb~ ro{z) ’ V—9tt\Y:g=^o{^)- (3.31) 

Therefore eq. (3.27) can be written in the simple form 


V-gttK\^^= 2 k 


(3.32) 


with constant k. 

That is, if in a given background spacetime we find a surface that satisfies (3.32), then 
we can ‘resolve’ this surface by replacing it and its interior with the static black hole with 
geometry (3.29), whose surface gravity is k. Let us remark that y/—gtt K in the limit to 
the horizon is, by definition, equal to 2k. The vector constraint is independent of p, but 
it is non-trivial that y/—gtt K takes the same value on the horizon and on . 

Using (3.13) we can write (3.32) more explicitly. For metrics on of the form (3.30), 
the square of (3.32) is 


772 (z) 


(^1 - ( D 77)2 + 



= 4k\ 


(3.33) 


When there is only one z coordinate, p = 1, this equation takes the form 

Uo2(z) / ^^^^2^ 


772 (z) 


(l - Y^TZ'izf + = 4^2 . 


(3.34) 


This is the version that we will employ in the examples in this paper. 

Finally, observe that in our construction of the black hole solution, K is positive when 
the radial normal points outwards of the horizon. This resolves any possible ambiguity 
about which side of in the background spacetime corresponds to the exterior of the 
black hole: if TLL > 0, the exterior of the black hole lies in the direction of the normal 
to 











3.3 Simple solutions 

We verify that the effective theory correctly reproduces known exact solutions. 


Schwarzschild-(A)dS black holes. In global AdS spacetime (extending to Minkowski 
and deSitter when < 0), 


ds — — ( 1 + ^ I dt + 


take a surface S at r = r(0), so that 

f^( 0 ) 




1 + v 


+ (dd^ + sin^ 6dQ.n+i) , 


(3.35) 


v,\e) = i + 


£2 ’ 


lee{0) = r{ey + 


(r-'(d ))2 

f(e )2 ’ 


1 + -^ 

It is immediate to see that a spherical surface 

r = rh 

solves (3.34) with the correct surface gravity 


1 / rf 


7^(d) =f(d)sin0. (3.36) 


(3.37) 


(3.38) 


The static planar and hyperbolic black holes in AdS are similarly easy to obtain. 

In a Minkowski background, l/£ = 0, we can also find the same solution starting from 

ds^ = —dt^ + dF + dF + Fdfl„+i . (3.39) 

Setting r = f{z), so that ^zz = 1 + {r'y and TZ = f, eq. (3.34) is solved by r = 
with K = l/(2r/j). Thus the Schwarzschild solution when 17 —>■ oo is obtained as a spherical 
soap bubble with mean curvature 

1 1 

(3.40) 


-K = —. 
n Th 


Black strings. Both in the Minkowski background (3.39) and in Poincare AdS, 

£2 

ds^ = (—dt^ + dF + Fdn„+i + dF) , 

the black string is just r = rh, with k = l/(2rh). 


(3.41) 


Intersecting black hole and deSitter horizons. Let us describe a less simple, lesser- 
known solution in the deSitter background,^ with = > 0, 

dr^ 


ds^ = - (1 - 77 ^ 2 ) df + 


1 - id 2^2 


+ (dd^ -|- siF 9dQn+i) 


(3.42) 


^We thank Jorge Santos for discussions that led to this example. 
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We change coordinates (r, 6) —>■ (f, 9) as 

sin 6 = Hr sin 9 , r cos 9 = r cos 9 , 


(3.43) 


so that (3.42) becomes 
ds^ = cos^ 9 (l — dt^ + — 
In this case, 


dH o \ 


+ sin^ 9d(l)" 


r = rh 

is a solution of (3.34) (where z = 9 and TZ{9) = r/^cos^), with 

This black hole is actually known exactly in all D [17]: it is the metric 


ds^ = cos^ 




with 


/ 


r 


1 


—fdt H-— + r d^n + 77^ ( d9 + sin ddcf) 


n—1 




- 


(3.44) 

(3.45) 

(3.46) 

(3.47) 

(3.48) 


which arises in the Kerr-deSitter family of black holes in the limit in which the equator 
of the rotating black hole touches the deSitter horizon and the configuration becomes 
static. The two horizons correspond to the cosmological and black hole horizons of the 
submetric (t,f, 0„), which extend separately along the coordinate 9 until they meet at 
the equator 9 = 7r/2. It is a simple exact instance of a horizon intersection. By the same 
reasoning as used in [2], in the large n limit the black hole of (3.47) becomes a ‘hole’ at 
f = = tq + 0{l/n) in the geometry (3.44). 

In the coordinates of (3.42), the solution (3.45) is the surface 


r = r{9) = 


Th 


— sin^ 0 (l — Td^r^) 


(3.49) 


If we write 


X = r cos 9 , y = r sin 9 , 


(3.50) 


then (3.49) becomes the ellipse 

^ + 77V = 1- (3.51) 

Since f(0) = < 1/H and f{7rl2) = l/H, the surface f{9) is an oblate ellipsoid that 

touches the deSitter horizon at its equator, in agreement with our interpretation above. 
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3.4 Effective stress tensor 

In order to match the near-zone solution to the far zone, they must share the metric (3.30) 
at their common boundary (in the overlap zone), and also their extrinsic curvatures. In 
this way the full geometry is smoothly glued between the two zones. This construction is 
equivalent to substituting the near-zone by a membrane with geometry (3.30) and with 
a stress-energy tensor given by the quasilocal stress-energy tensor of the near-horizon 
geometry measured in the overlap zone. Having this membrane stress-energy tensor as a 
source for the gravitational field in the background, ensures that the matching of the near- 
and far-zone geometries is C^. 

The overlap zone is the asymptotic region of the near-zone. Asymptotically at large 
p, and neglecting terms 0(6“^), the metric (3.29) becomes 

ds^ -VQ{z)df + rl{z)^ + (^jab{z) + ‘^r‘^iz)fab{z)^ dz'^dz’’ 

-\-TZ‘^{z) ^1 H-^ d^ln+i ■ (3.52) 

with the condition (3.27). When the scalar constraint (3.13) is satisfied, this is in fact a 
solution at all p: it corresponds to empty space (Minkowski or AdS) at large n. Therefore 
(3.13) is not specific to black holes. Instead, it pertains to the definition of near-zone 
asymptotics of large n gravity. 

The geometry (3.52) is then the reference metric required to define the effective stress 
tensor. In general there can be terms in the metric at large p of the type a^i,{z) p/n, but 
all the functions a^v{z) that appear are fully determined, up to gauge, by the Einstein 
equations in the overlap zone in terms of ro{z) and TZ{z) and their derivatives. In this 
sense, they are analogous to the hrst terms in the Fefferman-Graham (FG) expansion in 
AdS, which are fixed by the boundary metric. The terms at large p that behave like e“^, 
.., correspond to normalizable perturbations in the near-zone and are not determined 
by the boundary metric. The leading terms oc e~P will give the quasilocal stress-energy 
tensor. 

To leading order at large p, the quasilocal stress-energy tensor from (3.52) is the same 
for empty space and for the black hole. Thus this stress tensor cannot give any non-trivial 
gravitational effect on the background. Such effects come from the difference between the 
quasilocal stress energy tensors of empty space and of the black hole, i.e., 

87^GT^^, = -[K^,] + 6^,[K] (3.53) 

where the subtraction in = K^i, — is from empty space with the same 

asymptotic boundary metric. The subtraction removes the terms oc p at large p, and leaves 
those that decay like e~P (or faster). Typically, in the far zone we have e~^ ~ (ro/r)^, 
which is the fall-off of the gravitational field away from a localized source. 
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For the black hole solution (3.29), expanding at large p, neglecting terms 0{e and 


subtracting the background values we find 

277 

[K\] = —^e-P + 0{n^), (3.54) 

ro[z) 

[K\] = -2roiz)f\ e-P + 0{l/n ), (3.55) 

= + (3.56) 

[K] = 0{n°) . (3.57) 

It is easy to check that 

V^[KPa]-da[K] = 0{n^) (3.58) 


since this indeed follows from the vector constraint. It implies the conservation of the 
stress-energy tensor (3.53), which is given by 

277 

SttGT^ = - r^e~P + 0{n^) , 

mz) 

87rGr“b = C>(n°), (3.59) 

SvrGT^j = 0{vP), 

with the constraint that 


n(9bln7^)^^ - {dahin)5hT^j = r^(9,lnro) + 0(n°). (3.60) 

That is, at this order the values of T^a and T*j are not fully determined, since terms of 
0{'nP) in [K], which remain undetermined at this order, can make a contribution to these 
components. However, the combination in (3.60) does not suffer from this indeterminacy 
and this equation is required for the conservation of the stress-energy tensor at leading 
order. It is not surprising that the tensions (or pressures) play a subleading role: it was 
already observed in [2] that the pressure of black branes is suppressed relative to their 
energy density when D ^ 1. 

This effective stress tensor is enough to obtain the backreaction of the black hole on the 
background, by computing the linearized gravitational field created by this source. The 
specific values of T^b and are not important at this order, so one may choose them 
arbitrarily as long as (3.60) (which is necessary for a consistent coupling of the source to 
the gravitational field) is satisfied. While these components may be determined at the 
next order in the expansion, one way of choosing them is by requiring that, as is the case 
in all known black hole solutions, the ‘hoop stresses’ T^j vanish. Under this assumption 
we have 


SirCT^adblnTZiz) = ---5alnro( 2 :)e P + 0{l/n), 

mz) 


MGTj = 0 . 

In the case p = 1 we can solve (3.61) to find 

2 5^1nro(2;) 


8ttGT% = - 


ro{z) dzlnTl{z) 


e P + 0{l/n). 


(3.61) 

(3.62) 


(3.63) 
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3.4.1 Effective membrane source 


Since the stress tensor is homogeneous on we can integrate it over the angular 

directions. At large p this yields a factor Then the integrated energy 

density is 


{T\) = 


nO„+i 7^(z)"+l 
IhvrG rQ{z) 


+ 0{n^), 


(3.64) 


and the stresses are, as mentioned above, O(n^) and therefore subleading but must satisfy 


(^^)0bhl7^(z) = 


Qn+i 7Z(zr+^ 
IGttG ro(z) 


dalnro(z) + 0(l/n). 


(3.65) 


When p = 1 we can explicitly solve this to obtain 


(T%} = 


Qn+i'^iz^'^^djnroiz) 
IGttG ro(z) dzln7Z(z) 


(3.66) 


The energy density (3.64) admits a very simple physical interpretation; it is the energy 
density of a black brane of radius TZ{z), redshifted by the local redshift factor l/ro(z). 
The tension (3.66) is also the tension for that type of black brane, with an additional 
correction that accounts for the bending tension. 

The backreaction of the black hole on the background is obtained by taking a membrane 
with this stress-energy tensor, extended along the directions in the background, and at 
the origin of the 5”: by homogeneity in the sphere direction, in linearized gravity and at 
large n, we need not consider that the membrane extends on a sphere of finite radius but 
rather we can collapse it to zero radius. 

In appendix C we describe how the field created by this stress tensor correctly yields 
the large-n linearized Schwarzschild field, when it is interpreted as a solution of (3.34) in 
the background (3.39). 


4 Black droplets 

Now we construct solutions of (3.34) in the AdS spacetime (3.41) that are rather more 
complicated than in sec. 3.3. We seek surfaces that end at the boundary at z = 0 on 
a sphere 5"'”’'^, and extend a finite distance into the bulk at z > 0 until the 5'"'"’“^ shrinks 
to zero. They correspond to black holes localized at the AdS boundary. These were 
investigated in the context of cutoff-AdS/CFT (the holographic duality for the Randall- 
Sundrum-II braneworld) for the purpose of studying Hawking radiation of the dual CFT 
at planar (leading large N) order, and its backreaction on the black hole [18, 19]. In [13] 
they were reconsidered without the UV cutoff, in which case they are dual to the CFT 
in a fixed black hole background. After several years of controversy, the existence of such 
regular static black hole solutions seems to be settled after the numerical construction in 
[20]. Nevertheless, a simpler solution to the problem, such as afforded by the large D 
expansion, may be desirable. 
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We seek solutions for surfaces embedded in (3.41) in the form 


r = r(z). 


(4.1) 


Choosing a normal that near the boundary points outwards from in the direction of 
increasing r, a direct computation of the mean curvature to leading order in 1/n gives 



n 


z + rr' 


In this case eq. (3.32) becomes 



nz 


z + rr' 


= 2k . 


(4.2) 


(4.3) 


Note that i disappears from this equation, as it must, without having set it to 1. 

Eq. (4.3) is a non-linear equation that requires numerical integration, but it is conve¬ 
nient to first analyze its main properties. The equation simplifies for /I = 0 to 


rr' + z = 0, 


(4.4) 


which is solved by a circular profile, 

f‘^ + z‘^ = const. (4.5) 


The extrinsic curvature K of this solution vanishes: the curvature of the sphere is exactly 
cancelled by the cosmological curvature. Other simple solutions are: 

- the AdS black brane, with z = 1/(2 k) and f' —)■ -|-oo; 


- the black string solution, with r = l/(2«:). 


Eq. (4.3) is invariant under z —)> —z, and it requires that f' vanish at the boundary 
z = 0. Therefore, the surface always meets the boundary orthogonally. 

When extended into the bulk, we want to cap off smoothly at a finite value z = Zm 
where f = 0. The smoothness is guaranteed by eq. (4.3): for small f and large \f'\, it 
becomes 


rlr'l 


1 -|- 2kz 

which shows that if closes off away from the boundary, z > 0, it does it smoothly. 
The solutions of (4.3) are completely determined once the boundary radius 


(4.6) 


= r{z = 0) 


(4.7) 


is specified. Therefore, the solutions are parametrized by the surface gravity, i.e., by k, 
and by the radius at the boundary r^. These parameters can be varied independently, but 
only their product 

A = 2kr}j (4.8) 
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is invariant under the scaling symmetry of the background. Therefore we obtain a one- 
parameter family of inequivalent solutions labeled by A. 

Observe that, even if the near-horizon geometries for different values of A are all locally 
equivalent, they will match to different far-zone geometries. For instance, if we fix rj,, then 
as k varies we get boundary black hole geometries with the same horizon radius but with 
different radial dependence in gu- Ref. [13] proposed the existence of a one-parameter 
family of black droplets of this kind. Our construction realizes it in a natural way. 

When k is finite and non-zero we can rescale t,r,z such that we effectively set 


2 k = 1, 


(4.9) 


in which case the solutions are parametrized by A = r^. This choice is slightly convenient. 
With it, the two branches of solutions for f' in eq. (4.3) are 


z 1 ± + r^jl - z"^) 

r 1 — 


(4.10) 


Consider hrst the — sign solution. Substituting in (4.2) one finds that 

s^<-> = I ■ 

Given our choice of normal in (4.2), since K > 0 the criterion discussed at the end of 
sec. 3.2 tells us that the black hole exterior is in the direction of increasing r. Thus, the 
— sign in (4.10) gives black droplets. 

On the other hand, for the + sign in (4.10) we obtain® 

= (4.12) 

SO this gives solutions where the exterior is inverted relative to the previous case. Therefore 
these are not black droplets. It is less clear whether these solutions are relevant. We discuss 
them briefly in appendix D. 

Numerical integration of (4.10), with the — sign, yields the expected droplet-shaped 
solutions. In fig. 1 we plot them with fixed rj, = 1 and different values of k. 

We see that as k grows with fixed, the droplets extend more into the bulk. The 
parameter A ranges between 0 and 1, and as A —?• 0 the solutions approach (4.5). Eq. (4.4) 
is actually recovered from (4.3) in the limit in which both f and z go to zero at the same 
rate. At the opposite end, when A —)• 1, the solutions approach the black string. Near the 
black string limit, with A < 1, the solution at small z is close to a black string. 


2 kf = 1-(1 -A)e2 ^'^'. 


(4.13) 


Close to the cap it is instead described by (4.6). 

®To obtain this we assume that 0 < 2 < 1, which is appropriate since we want solutions connected to 
the boundary. 
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Figure 1: Black droplet solutions for fixed horizon radius at the boundary [vb = 1) and varying 
surface gravity k. The AdS boundary z = 0 is at the top, and z increases towards the bottom. 
The solutions correspond to A = 2krb = 0.01, 0.5, 0.95, 0.9995, with the elongation growing as A 
increases towards 1. 


It is not fully clear from this study which of these black droplets should correspond to 
the solution whose boundary geometry is the Schwarzschild black hole, although a plausible 
candidate is the solution (4.5). This, as well as the calculation of the holographic stress- 
energy tensor of these solutions, requires a more detailed analysis which we postpone to 
the future. 

Black droplets with black branes, and black funnels. In AdS one expects other 
classes of static black hole solutions that are related to droplets [13, 21]. Black droplets are 
often considered in conjunction with a black brane of infinite extent in directions parallel 
to the boundary. The two horizons are separated in the bulk, and in the large D limit the 
interaction between them is suppressed exponentiallly in D. Then we can obtain these 
configurations by simply superimposing a flat black brane and one of our black droplets. 
Their interaction through graviton exchange across the bulk is non-perturbative in l/H 
but in principle it is possible to incorporate such effects by computing the gravitational 
attraction between the two membranes with the effective stress-energy tensor of sec. 3.4.1. 

‘Black funnels’ can be regarded as the merger of a black brane and a black string (or 
a long enough black droplet) that hangs from boundary. Our effective equations do not 
seem to allow to obtain these solutions. The apparent reason is that the ‘shoulder’ at 
which the black string and the black brane are joined involves large gradients along the 
horizon, of order D, which fall outside the remit of the effective theory. 
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5 Effective theory: next-to-leading order 


At the next order we must take into account the p-dependence in K, in (2.8). It is convenient 
to separate it in the form 




(5.1) 


Here ro{z) is different than the leading-order function rQ(z) in (3.4), since it also contains 
1/n contributions. Since we will want to keep the horizon at /? = 0, we fix the ambiguity 
in this split by demanding that 

Jr(0,z) = 0. (5.2) 


We can obtain rg and 5r using eq. (A.5) in (2.8) and taking into account that, at this 
order, 


^ v 72 at DaroD“7?, 

— V A = n- 

N tqTZ 


Then we get 


(5.3) 


1 


(1 - (D7^)^) - + Lmr 

nl-i^ n l<, n tq I-c 


_ 1 2p + 3 f 1 
“o 4“ 


n 


1 


^2 


2 D'^n 3 DaroD“7^ 1 


n TZ n tqTZ 






(5.4) 


(recall that n and D are related by (2.9)) and 


(5r = ri(2:) ln(cosh(/9/2)), (5.5) 

with ^ 

niz) = || (l - rlr^DaTZDbn) , (5.6) 

and fab given in (3.20). In these equations Tl{z) and the metric ^abiz) with connection Dq 
are the ones obtained in the previous section, which already had absorbed 1/n contribu¬ 
tions from radial integrations. As with ro{z), they satisfy 

R{0, z) =n{z) , gabiO,z) = Jab{z) . (5.7) 

If we rescale p by a function f{z)+g{p, z)/n and shift 2; —)> z + h{p, z)/n‘^, we can reach 
a gauge with 

N{p, z) = ro{z). (5.8) 

We now write 

K = ———coth p + SK{p, z), (5-9) 

'^o[z) 

where the z-dependence ambiguity is fixed by demanding that 6K{0,z) be finite. The 
radial equation for 5K is 

dp5K + 2 coth p6K = —(5.10) 

ro 
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and its solution regular at p = 0 is 


5K = - 


ro sinh^ p Jq 
ri 


dp' sin]r{p')6r{p', z) 


ro 


where we define 


Next, setting 


Flip) = w 


sinh^ p Jq 


K\ = 


dp' In (cosh(/9'/2)) sinh^(/9'). 


D - 1 


+ 5K\ 


roiz) sinh p 

and using eq. (A.2) in (2.5) we obtain the equation 

o irr-t , ^1. ^I Flip) , 1 D„roD“7^ 

dp5K t + coth p6K t = — . , + ro ^-— 

ro sinh p \ Tq tq /•i 

which we integrate requiring regularity at p = 0, 


6KK,= 


ri F2ip) 

ro sinh p 


+ ro -n - 


1 DaroD“7^ 


ro7^ 


tanh(/9/2), 


with 


F2ip)= / dp'F^ip'). 


We can now integrate (2.6) to obtain gu = — V^(p, z). If we write it in the form 
V(/0, z) = Moiz) tanh(p/2) (l + -5Vip, z)^ , 


n 


with (5y(0, z) = 0, the result is 


where 


6V = riFsip) + 2rg ( ) lncosh(p/2) 




(5.11) 

(5.12) 

(5.13) 

(5.14) 

(5.15) 

(5.16) 

L 

(5.17) 

(5.18) 

(5.19) 


At this moment we can compute and 5K^j and use them in the vector constraint 
to obtain the NLO corrections to the effective equation. We have done this, but a quicker 
route to the equation is to impose directly the condition that the surface gravity on the 
horizon is constant, which must hold at all orders. Both conditions can be seen to be 
equivalent. 

Since we are fixing V(0,z) = \/oiz) and N(0, z) = roiz) at the horizon, the surface 
gravity is simply 

D-lVoiz) 


K = 


2 roiz) 


(5.20) 
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However, the effective theory uses functions of z in the asymptotic overlap zone, 1 <C p 'C 
n. In this region, 

Fs{p) = \ + 0{e-P), (5.21) 

so we find 

V(p, z) = Vo(z) (l + ^ J + 0{p/n) + 0{e-P). (5.22) 

In order to match this to a value 

(5.23) 

computed in the far-zone background, we set^ 

Vo(z) = I/(z) (l - . (5.24) 

Finally, in the asymptotic zone we have R(p, z)|^^= 7^(z) and ga6(p, z)|^^= '^ah{z)- 
Then /C(p, z)|^^= l/ro(z) and we can use its form (5.4) in (5.20). 


5.1 Effective equation at NLO 

Inserting (5.4) and (5.24) in (5.20) we obtain 


Ak? 


n{n + 1) 




2n 


2p + 3 
nP‘ 


2 _ 3 D^^roD^ ^ 

n IZ n roTZ rfi 


(5.25) 


where tq is the leading-order function given in (3.13), namely, 

r,>{z) = ,-„ . (5.26) 

Using (3.25) in (5.6), after a short calculation this equation can finally be written as 


Ak? 


z ^ + - 




1 1 /2p-k4 2D27^ 4DaroD“7^ 


n 




TZ 


tqTZ 




Da'^DfeT?.^ 

TV 


(5.27) 

which is the effective equation at NLO we sought. The last terms, with the Ricci tensor 
of 7af,, must be included in the second case of (2.14). 

In the particular case when p = 1 this effective equation is 


AV 




1 1 /6 






TIN = -2 + r ho - 27^^- (7^^)V - 47' 


,zz ' 0 


r'^TZ' 


with 


n \ V 


ro(z) = 


TZ 

TZ 


TZ 


tqTZ 


l_7^.(7^/)2+ 




(5.28) 


(5.29) 


^The subleading terms at 0{p/n) and 0{e~^) have to match between zones if they solve the equations 
of motion, since we have fixed the radial gauge and these terms do not correspond to any collective degrees 
of freedom. See the related remarks in sec. 3.4. 
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These NLO equations admit as solutions the ones discussed in sec. 3.3. Moreover, they 
also reproduce correctly the surface gravity of the corresponding exact black holes to this 
order. 

It is easy to see that if gu = 1 and the intrinsic curvature vanishes (or can be 

neglected), then ri is a constant (for any p) when the leading-order equation is satisfied. In 
this case, which typically includes surfaces in the flat Minkowski background, the effective 
equation at NLO is equivalent to the soap-hlm equation 

I constant (5.30) 

expanded to NLO in 1/n. In other cases, the covariant form of (5.27) is less straightfor¬ 
ward. 


6 Non-uniform black strings 


In the Minkowski background (3.39) the only solutions to the leading-order equation (3.34) 
are the spherical Schwarzschild black hole and the uniform black strings. However, at finite 
values of D it is known that non-uniform strings are possible [15, 16, 22]. Since they are 
absent from the leading-order large D theory, we seek them using the equation (5.30) at 
the next-to-leading order. Moreover, the wavelength of Gregory-Laflamme perturbations 
in the large D limit [1, 2] is ~ -y/n, which implies that non-uniform black strings must be 
sought among solutions of the second kind in (2.13). 

Following these remarks, we write the background geometry as 

ds^ = —dt^ H-h -|- r^dVln+i ■, (6.1) 

n 


and take 

7^(z) = 1 + . (6.2) 

n 

Here we have fixed the overall scale by setting the 5'"''’'^-area-radius of the uniform string 
to one. 

As a consequence of (6.2), we are limited to considering non-uniformities along the 
string of relative amplitude 0{l/n). Nevertheless, although small, these fluctuations are 
treated non-linearly. 

For (6.1) and (6.2), eq. (5.28) takes the form 


r”iz)+r{z) + r\zf = ro, 


(6.3) 


where Vq is related to the surface gravity k by 

= (6.4) 

Eq. (6.3) is non-linear in "P, but already its linear approximation at small V imme¬ 
diately reveals the existence of static, slightly non-uniform black strings with wavelength 
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Az ~ 2tt, which agrees with the value obtained in the small-amplitude perturbative anal¬ 
ysis of the Gregory-Laflamme problem at D —)• oo [1, 2]. 

Multiplying (6.3) by we obtain the first integral, 

V\zf = -V{z) +Vo + ^- . (6.5) 

It is useful to regard this equation as the classical mechanics of a particle in the 
potential 

V{V) = V + kie-^'^ , (6.6) 

with the coordinate z taking the role of time. Both positive and negative values of V 
are allowed as long as \'P\ <C n. When fci > 0 the potential is bounded below and has a 
minimum at "P = ^ ln(2fei). In this case there are solutions in which V oscillates between 
two turning points: these solutions are non-uniform black strings. When /ci < 0 the 
potential is unbounded below for negative P, which is not acceptable since then P will 
roll down to arbitrarily negative values, violating the condition |P| ^ re. Thus we only 
consider /ci > 0. 

We use the symmetry of eq. (6.5) under P —)> P -|- a, Pq —)• Pq + a, fei —)• to fix 

/ci = 1/2 so the equation is 

V'{zf = Po - V{z) + ^ (l - . (6.7) 

The constant Po corresponds to the energy of the particle and thus controls the range of 
oscillation of V{z) i.e., the amount of deformation. For Pq = 0 we find the uniform black 
string with V{z) = 0, and if Po ^ 0 we recover the small oscillations mentioned above. 

We are now interested in larger deformations, 1 <C Po re. The potential in (6.7) has 
two regions. To the right, where P > 0, it is dominated by the linear term ~ P. To the 
left, with P < 0, it is dominated by the exponential ~ In terms of particle motion, 

the particle will spend much more time {i.e., long extent in z) moving relatively slowly 
in the right side of the potential, than in the left side where it will move very quickly 
{i.e., very short extent in z) since the potential is very steep there. So the deformation of 
the string will have one large bulge (where P > 0) extending over a long distance, and a 
smaller neck (where P < 0) that extends over a small length. 

We can make this quantitative by solving (6.7) in an approximate way for large Pq. 
In the right side of the potential, for P > 0, we solve 

r'{zf ^ro + l-r{z) (6.8) 

i.e., a parabolic profile 

P~Po-|-- — -(^ — Zmax)"^ 5 Zmax = "^0 + ^ ' (6-9) 
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Figure 2: Profile of the non-uniformity of black strings for increasing values of the deformation 
parameter Vo = 1, 5, 10 (with Vo = Q corresponding to no deformation). Solid blue: numerical 
integration of (6.7). Dashed black: analytic approximation (6.9) [z > 0), (6.11) {z < 0). For 
clarity, we show only a half-period of the numerical results, and a full period of the analytical ones. 
The relative non-uniformity is MZ/TZ = 2V{z)ln and the proper length along the string is zj^fn. 


Here Zmax measures the distance from the point z = 0 of no-deformation, "P = 0, to the 
maximum deformation where V = Vmax ^ 'Po- In the left side, where P < 0, we solve 

r'izf ^Vo + ^{l- , ( 6 . 10 ) 

i.e., 

P{z) ~ 

where 

1 

Zmin — / = 

V^o + ^ 

is the length of the negative deformation from z = 0 to the minimum where P = Vmin ~ 
-^ln(2Po). 

The complete solution obtained by taking (6.9) for P > 0 and (6.11) for P < 0 ® gives 
an excellent approximation, as compared to numerical integration, to the entire profile of 
the deformation at large Pq. It also remains fairly good even for relatively small Pq (in 
particular if we keep the 1/2 offset for Pq in (6.9) and (6.11)). We present the two types 
of calculation in fig. 2. 

Let us now investigate the breakdown of the large-P expansion when the amplitude 

®The match at z = 0 is not completely smooth since the slopes differ at each side, but the difference is 
small for large Vo- 


arctanh a /1 — 


2Po + 1 2-v/Po 


ln(8Po) 


( 6 . 12 ) 


(2Po -|- 1) 1 — tanh (a / Pq ~\~ ~{z -\- Zmin)) 


( 6 . 11 ) 
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of the deformation grows large. Using the proper coordinate 


Z = 





(6.13) 


we find that, according to (6.9), the large bulge has extent 


(AZ)bulge ~ 2 



whereas from (6.11), the short neck has extent 


(6.14) 


(AZ)neck«^^^ln(8iPo). (6.15) 

If we consider a deformation at the limit of validity of our approximations, Vq ~ n, then 

( In Ti \ 

— \ . (6.16) 


Therefore, as the non-uniformity of the black string grows, it goes from having wavelength 
~ l/\/n and amplitude Vo/n = 0(l/n), to developing bulges of size Vo/n = 0(1), which 
resemble spherical black holes of 0(1) radius that are joined by thin necks. At these necks 
the approximation breaks down since the gradients are large, ~ n/lnn. 


7 Outlook 

In this article we have shown that an effective theory of black holes using membrane-like 
variables can be constructed in the large D limit. In spite of broad similarities to the 
membrane paradigm of [24], this theory differs from it in significant ways that may be 
worth spelling out. The basic idea in the membrane paradigm is to exploit large boosts 
near a null surface. Then it depends only the Rindler geometry. In contrast, the large-H 
effective membrane does not reside at the black hole horizon but at the outer boundary 
of the near-horizon region, which contains much more structure of the geometry than the 
Rindler limit. The large D theory is based on a neat separation between the black hole 
and the background, absent at finite D, which permits the explicit integration of short- 
wavelength degrees of freedom and gives an equation for the shape of the black hole. The 
membrane paradigm of [24], instead, does not yield any constraints nor information on 
the black hole shape. An equation for black hole horizons that is similar to (1.1) has been 
derived, in any D, in [23], but it is unclear to us if there is a connection to our approach. 

As in the case of the effective black hole theories of [9, 10], finding solutions using the 
large-H equations is considerably simpler than in the full Einstein theory. This is because: 
(i) the radial direction has already been integrated, resulting in a reduction by one in 
the cohomogeneity of the system, i.e., in the effective dimensionality of the equations to 
solve; (ii) there is only one degree of freedom (for static black holes), hence one equation 
to solve, instead of a system of coupled differential equations. 
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We have demontrated that two classes of non-trivial solutions are simply obtained with 
this method, but the detailed study of their physical properties has been outside the scope 
of this article. This includes the computation of the holographic stress-energy tensor for 
the black droplets, which requires to consider how the effective stress-energy tensor affects 
the far zone. The thermodynamic properties of the non-uniform black strings will be the 
subject of [25]. 

Clearly, there is much scope for extending the theory. The generalization to stationary 
black holes, which have an additional collective coordinate corresponding to a local boost 
along the horizon, will be presented elsewhere [26]. Including charge for the black holes 
is similarly straightforward [27]. The natural next step is the more general theory that 
includes not only spatial fluctuations but also time evolution on timescales At ^ vq/D — 
we understand that this is being pursued in [28] . 
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A Curvature on sections at constant p 

For a metric of the form 


gf.udx^^dx'' = -V2(p, z)dt^ + gab{p, z)dz'^dz^ + R^(p, z)dnn+i , (A.l) 


the Ricci tensor at any given value of p is (exactly in n) 



(A.2) 


(A.3) 


(A.4) 
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and the scalar curvature 


R = 


^(1-(DR)^)-2(„ + 1)(5^+. 


.ab 


9aR9bV\ 

RV ) 


-2^ + te)i?. (A.5) 


Here D and are the connection and Ricci tensor for the metric gab, and they are 

p-dependent. With our choice of fixing the z dependence at the horizon, up to NLO they 
coincide with the D and ^'^'^R^b of the metric 'jab both at the horizon and at the surface 
Tib in the asymptotic overlap zone. 


B Another derivation of the effective equation 


In this appendix we take a quicker route to deriving (admittedly, with some hindsight) 
the effective equation (3.33). Instead of solving the large D Einstein equations ab initio, 
as we have done in sec. 3, here we take an already known solution, from the large D 
limit of Schwarzschild(-AdS), and make its parameters become slowly-varying functions 
along directions parallel to the horizon. Then we impose Einstein’s equations to find 
the equations that the deformation functions must satisfy, and solve for any p-dependent 
correction to the metric required to obtain a solution. Although this approach is less 
systematic and complete than the one followed in sec. 3, it affords a quicker route to the 
derivation of the equations that can be useful in other instances. 

Consider the large D limit of a Schwarzschild black p-brane 

2 

ds'^ = —4k^rQ tanh^(p/2)dt ^-|—+ 5abdz°'dz^ + TZ^ (cosh(p/2))^dH„+i. (B.l) 


The solution is naturally parametrized by the sphere radius TZq, by the surface gravity, 
which we parametrize with the radius tq, and by the metric along the z directions. In 
order for this to be a solution, it must be that Tq ^ and that the metric along 

the z directions be flat. 

Now we promote these parameters to functions of z, taking the metric ansatz 




ds = -VQ{z)ta.nh^{p/2)dt + „ dp -h 7ab{z) + -6'jab{p,z) dz^-dz' 


n 


1 


n 


+TZ‘^{z){cosh{p/2))r^dQn+i 


(B.2) 


where we make the same assumptions about the dependence on p and z as in sec. 3. 
We have added a term 6^ab{p,z) expecting that it will be required in order to obtain a 
solution.® 

At this point we may proceed to impose the Einstein equations. However, we can 
anticipate that these equations will require uniform surface gravity on the horizon (for the 
Killing vector dt). So, in the interest of expediency, we fix this already by imposing (3.27). 


®Other terms could be added, but when solving the equations one easily sees that either they do not 
enter at leading order, or they must vanish. 
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For obtaining the Einstein equations we perform a Kaluza-Klein reduction on the 
large-dimension sphere 5”^^. Consider a metric of the warped-product form 


ds^ = g^y{x)dx^dx'' + 


(B.3) 


with /i, = 0,... ,p -|- 1 so the total dimension is p -|- n -|- 3. Hatted quantities are for the 

complete metric, unhatted ones for the (p -|- 2)-metric. The Einstein-AdS action is 


1 = j dP^^^^xy/^ (r + 


{p + n + 2)(j) + n + 1) 
£2 


j + n(n + 1 ) (^ ^-20 




e 

{p + n + 2){p + n + 1) 
£2 


(B.4) 


The Einstein equations along the directions x^ are (see app. A) 


G/ = G/ - (n + 

-6^'-' ^(n -0 l)(n -0 2) (d(p)^ + 2(n -0 !)□(/> - n(n + 


1 

+ 2 *^ 


_ (p + n + 2)(p + n + 1) ^ ^ 

- ' ^2 • 

The remaining equation is the one for the dilaton 


D(/) + (n + 1 ) (d(^)^ -ne ^ ^ 


(B.5) 


(B. 6 ) 


Now we apply this to the metric (B.2) and take the large n limit. To illustrate the 
method more clearly, we restrict ourselves to the case p = 1 . 

Erom (B.2) we have 

4> = InTl(z)-\ —ln(cosh(p/2)). (B-7) 

n 

It is straightforward to see how the different terms in the equations scale with re. To 
leading order one easily computes 


□(/. = 


re 


r2(p)cosh2(p/2) 
so that (B. 6 ) requires that 

1 




ro(z)2 7^(z)2 £2’ 


(B. 8 ) 


(B.9) 


i.e., the same as (3.13), which is indeed the effective equation that one solves in order to 
determine the embedding of the surface 


ds \^g= -4:k r^{z)dt +jzziz)dz + TZ^{z)dQn+i, 
in the background spacetime. 


(B.IO) 
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When this equation is satished, one easily sees that all the other equations (B.5) are 
also solved to leading order, except for the one for Ga^ = 0, since both da4>d^4> Vad^cj) 
are 0{n) (unlike all other off-diagonal entries) and thus enter at leading order. This 
equation requires that 

&lnr„ = a.lnR(l-^^|^) . (B.ll) 

When ro and TZ are not constant, this is solved by 

= ^^0 (J 2 - lncosh(p/2), (B.12) 


where in the last line we have used (B.9). We have also absorbed any integration function 
of z into jzz so that 6^zz{0, z) = 0. 

In this form, we have arrived at the metric (3.29) and the equation (3.33). 


C Far-zone Schwarzschild field from the effective stress ten¬ 
sor 


In this appendix we show how the effective stress tensor (3.64) yields the correct field 
for the large-n Schwarzschild solution, when this is obtained as a solution of the p = 1 
effective equation (3.39) of the form 


roiz) = Th , 


n{z) = rh 



(C.l) 


We want to recover the large n linearized field of the Schwarzschild black hole as the field 
created by a line source (a string) along the segment z G [—rh, rh], with an energy density 
(3.64) of the form 

/ ~ 2 \ {ri+l )/2 

(' - ?|) ■ 

At finite n such a string segment would give rise not just to a spherically symmetric 
monopolar field but also to higher multipoles. However we will see presently that when 
n —7- 00 these multipoles become subdominant. 

In the linearized approximation we need to consider the Newtonian potential created 
by this linear mass distribution. For a point mass source, the solution of 


= 167rGM(5(”+3)(x) 


(C.3) 


in the cylindrical coordinates of (3.39) is 

167rGM 1 

(n -b l)Q,n+2 (r^ -b z2)(n+l)/2 ‘ 


(C.4) 
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Hence, for the source (C.2) along the segment y G [—we have 


<h(r, z) 


nn+2''^J-r^ ^(^2 + (^_|)2)(r^+l)/2 • 


(C.5) 


We are in the far-zone outside and therefore at > rj^{l + 0{1 /n)). Therefore 

in the region near the tip of the segment along r = 0, we have \z\ — Vh > rhjn. due to the 
numerator, at large n the integrand is always strongly peaked around < 1/n and we 

can use the saddle point approximation. A simple way to do the calculation is by writing 


X 

n 


(C.6) 


so that 


« -lk±ir»rL r .fa 

’ Sl„+2 (r2 + (2_rfa/^)2)l”+'>''2 






1 


H„+2 ' \/n (^2 + ^2)(ri+l)/2 J_^ 

^n+1 


f 




r 


h 


(^2 ^2)(n+l)/2 


(C.7) 


where we have used that, at large n, 


f^n+2 ^ 



(C.8) 


Eq. (C.7) is the correct linearized field of the Schwarzschild black hole in cylindrical coordi¬ 
nates in Z) = n-|-4 dimensions. The essential point is that at large n the main contribution 
comes from the energy density near \z\ < Ijy/n, i.e., a point-like source, while the higher 
multipoles of field of the rest of the segment are strongly suppressed. 


D Black cavities? 

Taking the -|- sign in (4.10), the exterior of the black hole is reversed relative to black 
droplets. We refer to these putative solutions as ‘black cavities’, since they describe a 
spherically symmetric region of space enclosed inside a horizon. An example of what we 
mean by this is given by the deSitter foliation of AdS spacetime, 

This metric describes a region inside a horizon at p = £, which extends in the r-direction 
from the boundary at r —)> oo until it closes off in the bulk at r = 0. In (D.l), however, 

^°We thank Hong Liu for a discussion of this example. 
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Figure 3: ‘Black cavity’ (?) solutions for fixed surface gravity (k = 1/2). The AdS boundary 
z = 0 is at the top, and z increases towards the bottom. The exterior of the horizon is the region 
that extends towards smaller z (i.e., reversed relative to fig. 1), so at the boundary, the exterior 
of the black hole is the region 0 < r < r;,. The solutions correspond to ri, = 0.15, 0.5, 1, 2, 5. 


there are no large radial gradients ~ D close to the horizon, so this is generically different 
than the solutions we study here (even though p = t \s a, minimal surface with it' = 0). 

As in the case of black droplets, there is a one-parameter family of solutions to (4.10) 
with -|- sign, labeled by A in (4.8). They range from A = 0 for the half-spherical cap (4.5) 
to A —>■ oo for the black brane at z = l/(2k). In fig. 3 we plot some of them, keeping 
fixed the surface gravity and varying the boundary horizon radius r;,. As increases, the 
horizon grows in directions parallel to the boundary. The solutions resemble an AdS black 
brane with a horizon that is almost planar for r <C r^, but as r approaches r;,, it bends 
towards the boundary and meets it at r = r;,. 

It is not clear whether such solutions can be completed into far-zone geometries that 
are regular, in particular along the line r = 0 at the origin of the In principle 

one can fix arbitrarily the geometry at the AdS boundary, and it may be possible to write 
down boundary geometries with the required behavior, namely, a boundary geometry that 
is regular on and inside a horizon of finite radius r = rh, and which develops large radial 
gradients close to the horizon, as required for the large D expansion. If such a boundary 
geometry can be found, it must admit an extension into the bulk. If this is regular in the 
cavity enclosed by the bulk horizon, then one might legitimately characterize it as a ‘black 
cavity’. 
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